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Abstract 

The crawling motion of Dictyostelium discoideum on substrata involves a number of coordinated 
events including cell contractions and cell protrusions. The mechanical forces exerted on the sub- 
stratum during these contractions have recently been quantified using traction force experiments. 
Based on the results from these experiments, we present a biomechanical model of Dictyostelium 
discoideum motility with an emphasis on the adhesive properties of the cell-substratum contact. 
Our model assumes that the cell contracts at a constant rate and is bound to the substratum by 
adhesive bridges which are modeled as elastic springs. These bridges are established at a spatially 
uniform rate while detachment occurs at a spatially varying, load-dependent rate. Using Monte- 
Carlo simulations and assuming a rigid substratum, we find that the cell speed depends only weakly 
on the adhesive properties of the cell-substratum, in agreement with experimental data. Varying 
the parameters that control the adhesive and contractile properties of the cell we are able to make 
testable predictions. We also extend our model to include a flexible substrate and show that our 
model is able to produce substratum deformations and force patterns that are quantitatively and 
qualitatively in agreement with experimental data. 
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I. INTRODUCTION 



Cell movement over solid surfaces plays a key role in many every-day biological processes 
including embryogemsis, osteogenesis, wound heaung, and immune defense Q. For example, 
neutrophils chemotax towards a wound in order to prevent infection [2J . On the other hand, 
cell motility can play a significant role in disease; for instance, cancer cells spread out and 
intrude into healthy tissue by directed, active motion {3, 4, 5]. Hence, deeper insight into 
the biochemical and mechanical processes involved in cell crawling would be of great interest 
and importance. 

Despite their apparent differences, many eukaryotic cells share essential characteristics 
of their crawling motion 

a a. 

At the macroscopic level, cell motion often consists of 
several distinguishable phases: (i) extension of a membrane protrusion (pseudopod) at the 
leading edge, (ii) attachment of the pseudopod to the substratum, and (iii) detachment 
and subsequent retraction of the cell rear. These mechanical changes are mainly driven by 
polymerizing F-actin (protrusion) and myosin motors (retraction) 2fl- Both processes are 
regulated and synchronized in a spatio-temporal manner [8|. Additionally, in many hig her 



In 



organisms, detachment is regulated via biochemical changes of focal adhesions jil lloL 
other motile cells, on the other hand, focal adhesions are absent and a similar degradation 
mechanism has not yet been reported. 

Much of our understanding of cell motility has come from experiments on the social 
amoeba Dictyostelium discoideum which has been established as an experimental model 



system during the past decades 
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14|. These cells move rapidly (~ 10 /zm/min) and 



can be very sensitive to chemical cues. Also, the availability of a large variety of mutants 
allows quantitative insight into regulatory as well as mechanical aspects of cell motion. This 
paper is devoted to presenting a simple model for Dictyostelium crawling, with specific 
emphasis on the biomechanics of adhesive contacts between the cells and the substratum. 

One motivation for this study relates to recent force cytometry experiments in which 
the traction forces exerted by motile Dictyostelium cells chemotaxing on elastic substrata 
have been measured very precisely 15|, [l6|. The observed stresses range up to ~50Pa, 
giving rise to contractile pole forces, defined as the total force exerted in the front and 
back half of the cell, of ~90pN. Typically, the contractile forces are concentrated in spots 
of ~ fim size. These experiments also reveal a strong correlation between force generation 
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and morphological changes associated with the aforementioned three-stage cycle. Thus, the 
cell motion exhibits a mechanical cycle consisting of (i) a contraction phase, initiated by 
pseudopod attachment, in which the stresses increase; (ii) a retraction phase, in which the 
rear detaches and is brought forward. Consequently the cell shrinks and the stresses relax; 
(hi) a protrusion phase, in which the cell extends a pseudopodium in the direction of motion. 
At this stage, the pseudopodium does not exert noticeable forces on the substratum. The 
length of such a cycle is in the order of ~l-2min for wild- type (WT) Dictyostelium cells 
and ~4min in cells lacking myosin II, a motor protein responsible for cytoskeletal force 
generation 15]. The cell displacement of 15/zm per cycle is roughly constant. 

The exact nature of the adhesive forces between Dictyostelium cells and the substratum is 
not known. Most likely, the observed forces are transmitted through discrete contact foci on 
the ventral side of the cell. These foci are associated with F-actin rich regions which appear 
in spatial and temporal proximity to stress foci [ill, [l8]. Actin foci are spatially static but 
have a lifetime of ~20sec. Wild-type (WT) cells have ~5-10 foci. On the other hand, based 
on experimental results on cell detachment in shear flow, the number of microscopic adhesive 
bridges between cell and substratum is estimated to be ~10 5 [l9,[2(3]. Hence each adhesion 
focus is comprised of many bridges. 

It is reasonable to expect that to some extent, the cell speed should be controlled by the 
strength of attachment and the dynamics of detachment. Clearly, neither a non-adherent 
cell nor a cell that is unable to detach can move. However, between these extreme cases, the 
cell speed seems to depend only weakly on its adhesiveness [211] - In support of this, weakly 
adherent talin-null cells move with roughly the same speed as WT cells Mutants 
lacking myosin II move more slowly than wild-type cells, but cover the same distance per 
contraction cycle, i.e., the period of the cycle is increased. These cells do exhibit a much 
reduced motility on strongly adhesive substrata [2l|, as this combination places the cells 
in the extreme case of not having enough strength to contract against the adhesive forces. 
Finally, the over-expression of paxillin reduces the adhesion, but leaves the speed during 
folate chemotaxis relatively unchanged 221 ] . 

The importance of attachment / detachment dynamics for cell motility has been addressed 



in many theoretical studies 
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24 



25 



26] . These typically predict a strong dependence of 



cell speed on cell-substratum adhesiveness. Indeed, the prediction of an optimal adhesiveness 
is in excellent agreement with experimental findings on mammalian cells j^]. But, as just 
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discussed, the situation appears to be different in Dictyostelium. 

In these models, cell motion follows either from the protru ding activity at the front 
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26| or from asymmetric detachment during cell contraction 
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24 



. In the latter 



models, cell contraction is represented by internal forces acting on a visco-elastic cell body 
and the attachment/detachment dynamics are represented by an effective friction term with 



the substratum 



25, 



26 



28j] . However, the experimental observation of discrete binding sites 
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24j is more 



suggest that a representation by discrete, breakable springs as in Refs. 
appropriate. 

In this work we argue, that contraction takes place at a constant rate and that the cell 
speed is limited by the rate of detachment of the adhesive bridges. That is, the rate-limiting 
step in cell motility in this case is the peeling of the cell from the substratum. Based on 
stress patterns observed in Ref. [l5] we assume that cell detachment takes place mainly 
during the contraction phase and that protrusion forces contribute only a small amount to 
cell detachment. Therefore, our theoretical model of cell motion emphasizes the role of cell 
detachment during the contraction phase. Our model makes testable predictions about the 
cell speed under various experimental situations. These include crawling on substrata with 
varying adhesiveness and the variation of a number of cell-specific parameters. 



II. MODEL 



A. Components and Assumptions 

Our model focuses on the contraction phase of the motility cycle and does not explicitly 
treat the protrusive forward motion part of the cycle. Instead, as is shown schematically 
in Fig. [U the cell is assumed to maintain protrusive activity throughout the cycle during 
which cell material is constantly transported to the front. This notion is corroborated by the 
observation that over the entire cycle the cell speed shows only little variation (del Alamo, 
private communication). Then, we can define the cell speed as the displacement of the back 
of the cell at the end of the contraction phase divided by the cycle period. 

We assume that during the contraction phase, with duration r, myosins contract the 
cell body uniformly with a constant speed. This is motivated by (i) direct inspection of 
contracting Dictyostelium cells [15| and (ii) the observation that the in vitro myosin velocity 
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FIG. 1: Schematic cross-section of a crawling Dictyostelium cell illustrating the motility cycle. The 
part of the ventral surface that is in adhesive contact with the substratum is shown in red. (a) At 
the start of the contraction phase, the contact area is maximal, (b) During the contraction phase, 
the contact area shrinks while the cell is continuously transporting its body to the front, (c) At 
the end of the contraction cycle, the cell body is transported as far to the front as is allowed by the 
rear-most adhesions. Note, that it is assumed, that the protrusive force itself does not contribute 
significantly to the peeling of the rear, (d and e) During the relaxation phase, contraction stops 
and a full ventral adhesion area is re-established beneath the pseudopodium. 



is load- independent 29(. The assumption of a constant contraction rate is an essential 
difference to earlier work in which the cell is described as a one- dimensional network of 
contractile elements, each of which is exerting the same force on the nodes of the network 



23l |. Our choice is motivated by the fact that force balance implies that, when attached 
elastically to a substratum, the interior of such networks is largely stress-free. This is, 
however, in contrast to experimental observations which show that the stress field extends 
into the interior of the cell-substratum area, indicating that cells do not operate a contractile 
network with prescribed forces. We also assume that the cell contraction is not hindered 
by viscous stress of the surrounding medium. Indeed, as shown in Ref . [15] , the forces due 
to fluid drag on the moving cell are much smaller than the experimentally observed forces 



exerted on the substratum (~0.1pN vs. ~90pN, [151] ). Thus, the cell is always in a state of 
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FIG. 2: The cell-substratum contact area during different stages of the contraction cycle, (a) The 
start of the contraction cycle with the adhesion sites shown as red circles. The position of these 
sites is measured in a coordinate system with the center of the ellipsoid as the origin, (b) During 
the contraction cycle the cell contracts uniformly at a constant speed. The initial position of the 
adhesion sites is shown as open circles while the current position is indicated by a solid circle, (c) 
The end of the contraction cycle, with the remaining attached sites shown in green. 



mechanical equilibrium and the motion of the cell is quasi-static. 

Further, we assume that the cell is attached to the substratum via adhesive bridges. 
These bridges form with a fixed on-rate k + and dissociate with an off-rate which is both 
force and position dependent. The force dependence accounts for the fact that the potential 
barrier between bound and unbound state is lowered by an external force 2fj, l30|] . The 
)osition dependence incorporates a possible preferred detachment at the rear vs. the front 



coupler 

an 



23j |. These asymmetric adhesion properties are known to play a major role in mammalian 
cells, where focal adhesion complexes are cou ple d to intra-cellular pathways 31] . To our 
knowledge, and contrary to other systems 
been measured yet in Dictyostelium. 



32J, such a differential adhesion has not 



B. Rigid Substratum Model 

In our simulations, the adhesion area is represented by an ellipse with a fixed number 
(N) of randomly distributed sites that can adhere to the substratum. Their position Xi(t) 
at time t is measured with respect to the center of the ellipse. The amount of contraction is 
parametrized by the contraction rate A which can take on values between and 1 and which 
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is defined as A = (R — R T ) / R where R, R T are the semi-major axes of the ellipse at the onset 
and end of contraction, respectively. We divide the contraction cycle into 100 equal timesteps 
dt and at each timestep the new position of node % is given by Xi(t + dt) = (xi(t) — x m (t))(l — 
Xdt/r) + x m (t). Here, x m (t) is the location of the cell's center which is allowed to shift in 
order to ensure a vanishing net force on the cell (see below). The position dependence of the 
off-rate is chosen to depend on the component x along the direction of motion as follows: 

k.(x) = k^ b - [k., b - k.j]^—^- , (1) 

J x f -x b 

where Xf/b represent the front/back of the cell at the start of the contraction cycle and 
where are independent parameters of our model. The probability that a particular 

site adheres is given by the equilibrium value k+ +£ ^ ■ 

The attachments between cell and substratum are modeled by elastic springs with spring 
constant k s . In the case of a very rigid substratum we can ignore the deformations in the 
substratum. Then, the force on a single bond is given by Fi(t) = k s {xi(t) — aj°), where x® is 
the initial position of the bond. In principle, our prescribed displacement of the nodes can 
lead to a non-zero net force on the cell. To ensure a vanishing net force after each iteration 
we use the fact that the motion is quasi-static and allow the ellipse to shift and rotate. 
Specifically, we minimize the total energy of the springs at time t + dt, 

E s = irYl \R<p(xi{t) - x m {t)) + x m {t + dt) - x°f (2) 

where R v is the matrix describing a rotation by (p. An implementation of this minimization 
procedure revealed that the shift of the cell's center is small (< 5%) compared to the trans- 
lation of the cell for most of our model parameters and only became significant for small 
ratios of k-j/k-j,- To compute the resulting traction stress, a, we tile the substratum into 
0.05 i?x0.05 R squares and compute the total force per area for each tile. 



The force dependence of the off-rate is approximated by an exponential factor 



33|, 



k^xM) = exp (a^^) , (3) 

where we have defined the dimensionless parameter a = Rk s A/ {k b T). The molecular length 
scale A characterizes the width of the potential well which prevents the adhesive bridge from 
breaking and is of the order of 1 nm 33]. 
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Attachment of bridges to the substratum is assumed to occur with a force-independent 
rate constant k+. Binding rates decrease exponentially with the distance between membrane 

n 

and substratum [34| • Therefore we assume that attachment occurs only inside the contracted 
ellipse. We assume that k + is uniform across the contact area. The density of bridges on 
the membrane is assumed to be constant, such that the total number of available bridges 
that can attach at time t is proportional to the area of the contact area ~ N(l — Xt/r) 2 . 

The uniform contraction builds up stress and, consequently, a number of foci will detach 
during the contraction phase. To calculate the speed of the cell we first compute the smallest 
value of the x-component for all attached foci, x m i n (0), at the start of the contraction cycle. 
This corresponds to the left-most attachment point in Fig. [Th- At the end of one contraction 
cycle, we determine the focus with the smallest value of the x-component, x m j n (r) (left-most 
point in Fig. [lb). Then, the speed of the cell is given by (x min {r) — x min (0))/r. For each 
parameter set, we performed 1000 independent contraction cycles and parameter values were 
chosen such that there is at least 1 remaining attachment point. 



C. Elastic Substratum Model 



Traction force experiments that measure the position of fluorescent beads require the use 
of deformable substrata. The observed deformations are typically ~ 0.2/im comparable 
to the typical length of adhesion molecules 



35| . Under these conditions, the adhesive bridges 
cannot be treated as non-interacting springs. Rather, the elongation of a bridge under a 
prescribed cell contraction is influenced by the amount of substratum deformation caused 
by neighboring springs. 

To capture this effect, we simulated a deformable substratum with Young's modulus E 
as a two-dimensional triangular network of springs with spring constant k su b and rest length 
L. In these simulations, the initial conditions, the on and off rates of the cell nodes and the 
contraction procedure are the same as described above. Now, however, we need to compute 
the new positions of the triangular mesh vertices after each timestep. For this, we compute 
the total energy, given by 

h >. h Na 

m = -f EilvM - vM ~ L ? + i £(*«(*) - y^)f ■ (4) 

i,j a=l 

Here, yi(t) is the position of the i-th triangular mesh vertex at time t. The first sum in 



Eq. H] extends over all pairs of neighbors in the triangular grid and the second sum runs 
over the substratum nodes that are coupled to N a adhesive springs. For simplicity, we 
have chosen boundary conditions in which the position of the substratum boundaries are 
fixed. Minimization of Eq. [4] directly yields the new positions of the vertices and, thus, the 
deformation pattern of the substratum. The force exerted on the attachment point yi{t) by 
the cell node Xj(t) can be calculated as F^(yi, i) = k s (xj(t) —yi(t)). The total force F(yt, t) 
on each attachment point is then the sum over all nodes j connected to this point. These 
point forces are related to the local applied stress via 

0vz{Vi,t) = -J^-JY~ ' v = x,V- ( 5 ) 

Note that our choice for the boundary condition will lead to non-zero net forces on the 
cell. We found, however, that for a substratum of sufficient size (4Rx4R) the net force is less 
then 5% of the pole force. Of course, by repositioning the cell after each time step we could 
ensure a vanishing net force even in the case of fixed boundaries. Furthermore, choosing 
periodic boundary conditions for the substratum will also guarantee a vanishing net force 
on the cell. We found that the resulting force pattern differs only slightly from the force 
pattern generated using fixed boundaries, demonstrating that the results are insensitive to 
the precise details of the numerical algorithm. 



D. Parameter Estimates 



Throughout the paper we will use a default set of parameters that were obtained, where 
possible, from experimental data. The shape of the cell is characterized by a long semi-axis, 
taken to be R =10um, and an aspect ratio 1:4. Based on movies shown as supplemental 



material to Ref. [15j and direct measurements of the adhesion area in Refs. [36|, [37|], we 
assume that the (WT) cell contracts by 50% of its length, corresponding to A = 0.5 in our 
simulations, during a contraction period of r=l min. 

For the number of adhesive bridges we followed Refs. 34j and chose N = 200. Note, how- 
ever, that our results do not depend on N as long as we rescale the other model parameters 
appropriately. Specifically, if N — > fiN we need to rescale k s and A as follows: k s — > k s j\i 
and A — > /zA. The off-rates are estimated in models of shear flow induced detachment 



20. 



34j and at the back we take = 1 xlO 2 /sec. As discussed before, there is no clear 
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data on the possible maturation of adhesion sites in Dictyostelium and we have arbitrarily 
chosen the off-rates at the front to be equal to O.bk-^- The force dependence of the off-rate 
in Eq. [3] is determined by the dimensionless parameter a which we have chosen to be 125. 
This parameter is a combination of the rupture width A of the molecular bond and the 
adhesive spring constant k s . We have chosen the latter to be k s = 1 x 10~ 4 N/m, which is in 



the range of experimental values 35j, and A ~ 0.5nm 33]. Finally, the spring constant of 
the deformable substratum was estimated using the experiments results in Ref. [15| . There, 
the pole force was found to be F p ~ 200 pN while the deformation was u ~ 0.2/zm, leading 
to k sub = F p /u = 1 x 10~ 3 N/m. 



III. NUMERICAL ANALYSIS AND RESULTS 



A. Rigid substratum 

With the above choice of parameters, we performed 1000 contraction cycle simulations. 
At time intervals dt = O.Olr the distribution of displacements 7, i = 1...N was 

stored. The displacements u x and u y are directly related to the traction forces exerted on the 
substratum via Fj = k s Ui, i — x,y. Fig. [3] shows the time evolution of the stress averaged 
over 1000 individual runs. Here, forces were summed up in bins of size 0.05 i?x0.05 R. 
Note, that the ellipses in our simulations correspond to the adhesion area which does not 

nn 

necessarily correspond to the experimentally determined cell outline |15l . 1161 ] . 

Fig. H] shows the force distribution averaged over time for different choices of model 
parameters. Averaging was done by scaling individual time frames such that the contracted 
ellipses fall on top of each other. The top pattern corresponds to the average stress pattern 
for the default parameters. In each row of images, we have varied one of these parameters 
and have plotted the stress pattern using a gray-scale with black corresponding to large 
stresses. Surprisingly, variations in the amount of contraction A and the relative adhesiveness 
k_j/k_^ have only little influence on the stress pattern. Rather, it depends more strongly 
on the basal detachment/attachment dynamics via k_^ and k + , along with the molecular 
length scale A and the spring constant k s . 

In Fig. [5]we plot the dependence of the pole forces on the model parameters as a function 
of time during a contraction cycle. The pole force at the back, Ft,, is defined as the total 



10 




t=0.5 
t=0.6 

?=o, 

1=0.8 

1=0.4 " -1=0.9 

FIG. 3: Average traction stress patterns over 1000 simulation runs with time expressed in units 
of the contraction cycle. For the purpose of averaging the distribution maps were tiled into 0.05 
i?x0.05 R squares. The stress is shown in a gray-scale with black corresponding to a traction stress 
of a ~4.4fc s i? _1 . At the beginning of the contraction cycle (t=0) no force is exerted. The outer 
ellipse indicates the original position of the cell and the inner ellipse indicates the current adhesion 
area. 

force exerted in the direction of motion, i.e. 

F b = k s J2 « • (6) 

u x >0 

Similarly, the pole force at the front, Ft, comprises all forces which point into the negative 



x-direction. Our definition of the pole forces differs from the one in Ref. 15j], where pole 
forces are defined as the overall forces transmitted at the attachment regions in the front 
and back halves of the cell. In each graph we used the default parameter set and varied one 
value as indicated in the legend. 

In Fig. O we compare the dependence of the cell speed on four model parameters. In 
our model, this speed is determined by the amount of retraction at the rear of the cell per 
contraction cycle: during the protrusion phase the ellipse representing the protruded cell 
outline is moved such that the rear coincides with the last remaining attached focus. The 
actual forward motion is accomplished throughout the contraction and protrusion phase 
(see Fig. [1]). Note that even for a symmetric detachment the cell can move forward. Again, 
we varied one parameter value with the remaining parameters fixed at the default values. 
Finally, in Fig. [7] we plot the average pole force during a single motility cycle as a function 
of the off-rate As expected, the pole force decreases as the off-rate increases. 
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k f =0.1k_ b k ,f=0.9k 




A=0.1 A=0.9 




FIG. 4: The traction stress a, averaged over an entire contraction cycle, for different sets of 
model parameters. The stress is plotted using a logarithmic gray scale with black corresponding 
to \a\ ~ Q.bkgR^ 1 and white corresponding to values |<t| < 6.5 • lO^kgR^ 1 . The time averaging 
was achieved by rescaling and overlaying the contracted ellipses. The upper pattern corresponds 
to the default set of parameters: /c_ j 6=6-10 _1 t _1 , A;_ j=0.5k^^, k + =6-10~ 1 r~ 1 , a=125, A = 0.5, 
and N = 200. For this set of parameters, the maximal stress is ~0.6fc s ii _1 . In each row one model 
parameter is varied while keeping the remaining parameters fixed. 

B. Elastic Substratum 

Fig. shows a time series of the displacement pattern for an elastic substratum. The 
displacement is shown using the indicated color scale and the computed maximal displace- 
ments ~ 0.02i? ~ 0.2/im are in good agreement with experimental results. Due to rapid 
detachment, substratum deformations vanish shortly after onset of contraction. Fig. [HJo 
shows the corresponding stress distribution [a 2 xz + cr^) 1 / 2 which shows two distinct peaks 
at the front and the back corresponding to regions of maximal displacement in Fig. [HK 
The computed maximal stresses (~ 4k s R~ l ~ 40Pa) are similar to the ones observed in 



12 




0.2 0.4 

time [t] 



0.2 0.4 0.6 0.8 

time [t] 



FIG. 5: Average pole-forces as function of time during one contraction cycle. The default parameter 
set is used and the parameter value indicated in the legend is varied. 




FIG. 6: The dependence of the cell speed on one out of the six model parameters is shown. The 
remaining parameters are fixed at their default values. 



experiments (~ 50Pa 15]). 




k-,b K 1 ] 



FIG. 7: The average pole force exerted during one contraction cycle as a function of 5. 
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FIG. 8: Time evolution of the displacement pattern (a) and stress pattern (b) of a deformable 
substratum. Displacements are given in units of R using the displayed color scale while stresses 
are given in k s /R, as shown in the color scale. In this simulation, the default parameter values 
were used together with an effective spring constant k su b = Wk s = 1 x 10~ 3 iV/m. 

IV. DISCUSSION 



In this paper we have presented a mathematical model for cell motility motivated by 
experimental observations of the motion of Dictyostelium cells. The emphasis of our model 
is on the interaction between the cell membrane and the substratum on which the cell is 
crawling while the actual cell deformation and translation are not explicitly taken into ac- 
count. There are several distinct differences between our approach and previous modeling 
studies. The studies carried out by Lauffenburger and co-workers 23|], for example, consid- 
ered a one-dimensional cells with only a handful of attachment points. These points were 
connected through springs that are exhibiting a prescribed force. In our model, on the other 
hand, the foci are moving with a constant contraction rate. This choice was motivated in 
part by the observed stress and force patterns in traction force experiments. These experi- 
ments demonstrate that the forces are maximal within the contact area. In a model where 
the inter-foci springs exert a fixed force, the force field within the interior of the contact area 
will be very small and concentrated at its boundary. Furthermore, experiments on TalinA- 
cells 15| demonstrate that cells with a vastly reduced adhesion move with roughly the same 



velocity as wild-type cells. A prescribed force model would predict a strong dependence of 
the cell's speed on the adhesion strength. 

Another major difference is the presence of dashpots, representing the viscous nature of 
the cell's cytoplasm, in the earlier models. These dashpots play an important role if one 
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prescribes the force exerted by contractile elements. Here, however, we prescribe the con- 
traction velocity which alleviates the need for an explicit modeling of the viscous cytoplasm. 
The estimation of the contraction speed we used in our simulations, ~ 10/im/min, is based 
on direct experimental observations. However, typical in vitro myosin velocities measured 
in motility assays are ~ 10 — 20 times higher than the experimentally observed cell speeds 



29l . [38|. The in vivo velocity is not known but will likely be of the order of the contrac- 



tion speed. The mechanism responsible for this significant slow down is, to our knowledge, 
unclear. One possibility is that in vivo the disordered structure of the actin-myosin cortex 
hinders a rapid contraction. Also, the viscosity of the cytoplasm may play an important 
role in limiting the myosin contraction speed. A final difference is that our two-dimensional 
model explicitly takes into account the displacement of the rear. 

Another class of models describe the cell as a gel, with visco-elastic properties 
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26 



28]. 



Contrary to our model, these studies prescribe the protrusion of the cell and do not focus 
on the contraction mechanism. In these models, the adhesion has a front-to-back gradient 
and is represented by an effective friction force. Thus, they are unable to address the role 
of contraction on the detachment of crawling cells. 

Most of our results are obtained assuming that the substratum is rigid, corresponding to 
a typical experimental set-up where cells are crawling on glass surfaces. In this case, the 
displacement of the adhesion proteins is much larger than the displacement of the attach- 
ment point at the substratum. Thus, the force field exerted on the substratum is simply 
determined by the forces on the adhesion proteins. As expected, this average traction force 
varies during the contraction cycle and reaches its maximum shortly after the start of the 
cycle (Fig. [3]). The pattern observed in Figs. [3] and H] can be explained by realizing that 
in our model stress is generated by a prescribed isotropic contraction. This leads to radial 
increase of stress at the adhesions which, in the absence of binding/unbinding dynamics, is 
given by the geometry of the contraction only. Thus, in our model the binding sites at the 
center of the adhesion zone are always almost stress free, resulting in the observed pattern. 

The default set of parameters of our model were based, where possible, on experimental 
values. To examine the effect of these parameters on the force patterns, we have systemati- 
cally changed one while keeping the remaining parameters fixed (Fig. H]). The stress pattern 
depends strongly on the molecular length scale A with the stress increasing for smaller val- 
ues of A. This parameter determines the off-rate of the bridges (Eq. [3]) and for small values 
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of A, this rate becomes small. Correspondingly, the force per focus becomes large, leading 
to large stresses shown in Fig. |H Note that the parameter k s also controls the off-rate. A 
change in k s , however, does not change the force pattern as dramatically as a change in A 
since this parameter determines the force per bridge as well. 

The relative adhesiveness k-j/k-^ measures the asymmetry in the adhesion strength 
between the front and the back of the cell. Such an asymmetry is essential to the motility 
of mammalian cell but its role in Dictyostelium movement is unclear. Variations in the 
amount of the relative adhesiveness k-j/k-^ have only little influence on the magnitude of 
the observed stress pattern. The pattern, however, becomes more asymmetric as k-j/k-fi 
decreases. Clearly, a larger off-rate at the front than at the back will lead to a higher 
concentration of attached bridges at the front and thus a larger stress in the front half of 
the cell. 

The parameter A describes the amount of contraction. In the absence of detachment, a 
larger contraction would lead to an increase in the elongation of the bridges and a larger force 
per area. However, the increased force on the foci will lead to an increase in the detachment 
and, as can be seen from Fig. HI these two effects compensate and lead to a slightly smaller 
time averaged stress for larger contractions. 

The on-rate k + describes the re-attachments of foci and increasing the value of k + will 
result in an increase in the number of attached foci during the contraction cycle. Thus, the 
force per area increases for increasing values of k + , as is evident from Fig. HI The off-rate 
on the other hand, determines the detachment dynamics of the foci. A higher value of 
leads to a smaller number of attached foci and thus a smaller force per area. 

The pole forces, defined as the sum of all the forces parallel or anti-parallel to the direction 
of the motion, increase rapidly and linearly at the start of the contraction cycle (see Fig. 
[5j). This linear behavior can be understood by realizing that during the initial contraction 
period, the force dependence of the off-rates is insignificant and the number of bridges stays 
roughly constant. Since the force on each adhesion is proportional to the contraction ratio, 
the pole force increases linearly. Once force induced detachment becomes significant the 
bridges begin to break and the pole force starts to decrease. The maximum pole force, and 
the time at which this maximum is reached, depend on the model parameters (Fig. [5]). In 
particular, the maximum value increases for smaller values of k-j, (Fig. [5^). After all, small 
values of the off rate lead to larger displacements and, thus, larger forces. Furthermore, the 
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pole force increases for larger values A (Fig. [5b) which can be understood by realizing that 
small contractions lead to small displacements and thus smaller pole forces. 

Using our model, we are able to vary systemically each parameter and determine the de- 
pendence of the speed on this parameter. The results (Fig. \5§ can be viewed as experimental 
predictions even though we realize it might be difficult to vary some of these parameters in 
experiments. In particular, it is not always obvious which adhesion parameter is probed in a 
certain experiment and how the parameters are changed in a certain mutation. For example, 
the reduced adhesiveness of TalinA- mutant may result from an increased off-rate or from a 
smaller total number of adhesive sites. Surprisingly, we find that the speed is only weakly 
dependent on the relative adhesiveness k_j/k_^. This is in contrast to previous models 
where the speed depends critically on this ratio. Our model assumes that the protrusion 
is decoupled from the contraction cycle (Fig. [1]). Thus, our speed is mainly determined by 
the peeling velocity of the back and can be significant even for uniform off rates. Note that 
for small relative adhesiveness it becomes important to ensure a vanishing net force through 
a re-orientation of the cell outline. Without this re-orientation the cell's speed would be 
purely given by the off-rate at the back and would be constant for all values of the relative 
adhesiveness. 

As expected, we find that the cell speed increases for increasing values of the contraction 
rate A (Fig. [6b). After all, in the limit of vanishing contraction rate the speed approaches 
zero while for maximal contraction rate the speed reaches a maximum. Furthermore, we 
find that high on-rates decrease the speed (Fig. [6b). For high values of k + , adhesive bridges 
are deposited at rates that are higher than the detachment rates, limiting the cell's speed. 

Contrary to previous studies, we find that the speed does not depend strongly on the 
off-rate k-^ (Fig. [6H). Of course, the speed will approach for very small values of this 
off-rate where the foci will remain attached to the substratum. In this limit, we expect 
that our constant contraction speed assumption is no longer valid and that the forces on 
the myosin motors are large enough to lead to stalling. For large values of the off-rate, 
all foci will detach and we have only considered the range of values for which at least one 
focus remains attached. In fact, in this limit the weakly adherent cells can exert only small 
forces on the substratum, see Fig. [7J Hence, for sufficiently large the traction force that 
balances the viscous drag of the protruding cell (~ O.lpN [15|) exceeds the detachment force. 
For approximately symmetric cells, force balance then implies that a forward protrusion is 
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accompanied by a backward motion of the same order. Hence, there is no net motion for 
sufficiently large k-j,. For the parameter range studied, the traction force is always sufficient 
to support protrusive forward motion, see Fig. [7J 

Our finding that the cell speed is roughly constant for a large range of values of adhe- 
sive forces is in agreement with recent experiments in which the stress patterns of crawling 
Dictyostelium cells were examined. These experiments show that the cell motion can be 
described by a contraction-relaxation-protrusion cycle. Thus, the cell's speed is determined 
by the ratio of the displacement per cycle and the period of this cycle. TalinA- cells exhibit 
a drastically reduced cell-substratum adhesion but were found to have the same cell speed 
as wild-type cells, with an identical period and, thus, identical displacement. Of course, 
two data points cannot rule out a significant dependence of the cell speed on the adhesion 
strength and a definitive test of our model would be to examine the cell speed for different 
mutants. One candidate would be cells in which the expression level of PaxB, the Dic- 
tyostelium orthologue of paxillin, is altered. Both PaxB- cells [39|] and cells in which PaxB is 
overexpressed 22] exhibit a decrease in cell-substratum adhesion. The cell speed in cAMP 
gradients is reduced in PaxB overexpressed cells and is increased in PaxB- cells. Interest- 
ingly, the cell speed in folate gradients is largely independent of the expression level of PaxB 



221 ] and D. Brazill, personal communication). This might indicate a PaxB role in the pe- 



riodicity of the motion cycle, which would effect the cell's speed. A more detailed analysis 
of these mutants that can measure force patterns and motility cycles would be interesting. 

A quantitative comparison with the experimentally obtained stress patterns is only pos- 
sible if we take into account a deformable substratum. After all, these experiments measure 
the displacement of fluorescent beads embedded in the substratum and require significant 
movement of these beads. Thus, our model assumption that the displacement of the sub- 
stratum is negligible compared to the stretching of the adhesive bonds is no longer valid. To 
compare to experiments, we have extended our model and have explicitly simulated a trian- 
gular spring network, representing the substratum. This extension renders the simulations 
computationally more demanding and we have only performed a limited set of simulations 
(Fig. [8]). Using experimental values characterizing the substratum, we found that our results 
show a quantitative and qualitative agreement with the experimentally observed stress and 
strain patterns. For our experimentally based parameter values we obtained a maximum 
displacement that was comparable to the one observed in experiments (~ 0.2/im). Further- 
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more, the computed peak stress is similar to the experimental peak stress: ~40 Pa vs. 
-50 Pa. 

In summary, we have presented a simple model for the motion of Dictyostelium cells. We 
have shown that this model can produce a number of experimentally verifiable predictions 
and can be extended to include deformable substrata. Our strongest prediction, that the 
cell speed is largely independent of the value of the adhesive forces, should be testable using 
force cytometry experiments. Our model focused on the cell-substratum interaction and 
ignored the protrusion phase of the motility cycle. Extensions that include intra-cellular 
signaling pathways that drive cell deformations are currently under investigation. 
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